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PROCEDE DE CALCUL DE REALISATIONS MAILLEES D'UN 

RESERVOIR 

La presente invention conceme le domaine de la g6ologie et plus 
5 particulteremen t de F exploration p etrol i dre. 

II est connu, notamment dans l'exploration p6troli6re, de determiner la position 
des reservoirs petroliers a partir des resultats de mesures geophysiques effectuees 
depuis la surface ou dans des puits de forage. Ces mesures impliquent typiquement 
remission dans le sous-sol d'une onde et la mesure des diverses reflexions de Tonde 

10 sur les structures geologiques recherchees- surfaces separant des mat&iaux distincts, 
failles, etc (technique de la sismique reflexion). D'autres mesures sont effectu6es 
depuis des puits. On envoie alors dans le sous-sol 5 des ondes acoustiques, des 
rayonnements .^arnma on des signaux 61ectriques. Ces techniques impliquent le 
traitement des mesures pour reconstituer un modele du sous-sol. Le traitement des 

1 5 mesures est aff<;ct6 d'incertitudes de diverses natures* Pour la sismique reflexion : 

Une premiere incertitude est l'incertitude de points. On appelle points pour des 
mesures geophysiques le suivi sur les images ou sur les sections sismiques 
successives fournies par traitement des mesures d'un marqueur sismique, par 
exemple un horizon geologique, pour definir une surface. II est possible pour 

20 l'op£rateur de confondre des marqueurs voisins ou encore d'associer des marqueurs 
qui ne font pas partie d'une mSme surface. L'incertitude de points peut etre de 
quelques m&tres k quelques dizaines de metres; la valeur de quelques metres 
correspond peu ou prou a la resolution, i.e. k la largeur des impulsions des signaux 
€mis dans le sous-sol; la valeur de quelques dizaines de metres correspond k une 

25 erreur de positionnement d'un marqueur sismique. 

Une autre incertitude est l'incertitude sur le champ de vitesses dans le sous-sol. 
La remise en profondeur des mesures sismiques effectives dans le domaine temps 
(cette remise en profondeur peraiet de passer d'une representation temporelle k une 
representation spatiale; la representation temporelle est celle foumie par les 

30 instruments de mesure - qui mesurent typiquement les reflexions de Tonde emise. La 
representation spatiale est obtenue k partir d'hypothfeses sur le champ de vitesses dans 
le sous-sol. Ellc suppose en effet des hypotheses sur la nature des materiaux traverses 
et sur la vitesse de propagation de I'onde dans ces mat&iaux. Or une hypothese de 
champ de vitesse est affectee d'une incertitude. Cette incertitude sur le champ de 

35 vitesses peut conduire k des incertitudes de position de Tordre de 50 mfetres. 

Une autre incertitude est Tincertitude due a la migration. On appelle 
"migration" dans ce domaine Poperation qui permet de repositionner correct ement les 
6v6nements dajis Fespace. L'incertitude sur la migration est fonction du contexte 
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g6o1ogique; elle peut etre quasiment nulle pour des structures tabulaires ou quasi- 
tabulaires et atteindre des centaines de metres pour des structures complexes, comme 
des domes de sel avec des parois verticales ou quasi -verticales. 

La presence de structures particulieres - des failles par exemple - attenuant la 
5 reflexion ou provoquant des diffractions de 1'onde emise provoque encore d'autres 
incertitudes. 

n est possible de relever grace aux mesures effectuees dans les puits les limites 
de structures geologiques. L'information ainsi obtenue est aussi entachee 
d'incertitudes; oes incertitudes proviennent essentiellement des mesures de position 
10 du puits. 

Des incertitudes de calage sont dues au report des limites de structures 
g6ologiques me$ur6es dans le puits sur toute section sismique intersectant le puits; 
elles proviennent pour Fessentiel des problemes de resolution. 

FR-A-2 797 329 propose une methode d'eiaboration de cartes de probability de 

15 position d'impact d'un puits. Ce document contient une discussion des incertitudes 
liees aux migrations temps et profondeur et des erreurs de position correspondantes. 

Compte tenu de ces incertitudes, le traitement de mesures g6ophysiques peut 
conduire a proposer plusieurs realisations des structures du sous-sol. Pour des 
reservoirs — qui sont les structures recherch6es par Pexploration petroliere - on 

20 propose des surfaces representatives du toit du reservoir, de la base du reservoir et 
d'interfaces entre des couches intermediaires. Une realisation est alors constitute d f un 
ensemble possible de surfaces representant le reservoir, correspondant a un ensemble 
d*hypotheses. Typiquement, en exploration petroliere, on fburnit jusqu'a 300 
realisations correspondant aux differentes hypotheses possibles. A titre indicatif, la 

25 fourniture d'une realisation (dite realisation de reference) a partir des resultats de 
mesures geopbysiques est une operation qui mdme avec l'assistance d'un outil 
logiciel, corres])ond a un travail considerable, pouvant atteindre 1 mois/homme pour 
un specialiste routine; lorsque la sismique a ete interpret6e, la construction d'une 
description surfacique d'un reservoir peut prendre une semaine ou plus. La collection 

30 des incertitude*, represente egalement une semaine, leur integration dans un logiciel 
une journee et 1e calcul des realisations quelques heures. . Sur une superficie globale 
de 20 km 2 , une surface peut etre representee par des triangles de 1 00 m de cote, soit 
done quelques milliers de points pax surface. Une realisation peut impliquer outre les 
surfaces du toit et de la base du reservoir, 1 & 5 surfaces intermediaires. 

35 FR-A-2 652 ISO d6crit un proc6d6 de modelisation de surfaces. II est suggere 

d'utiliser une description informatique maillee des surfaces. A certains des nceuds de 
Ja surface maillee sont associees des contraintes — representatives par exemple de la 
position connue de certains ncEuds ou d'une incertitude donnee sur la position de 
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certains nceuds. La description infomiatique de la surface est ensuite adaptee aux 
contraintes, en utilisant une methode appeiee "Discrete Smooth Interpolation" (DSI). 

Par ailleurs, pour representer les proprietes petrophysiques d'un reservoir, on 
utilise une grille ou module voluxnique maille, comnie decrit dans le brevet US-A-4 
5 821 164. Le procede d6crit dans ce document comprend des etapes de construction 
d'un plan separe pour chacune des surfaces critiques limitant le volume et de 
construction de couches de cellules en fonction de la structure stratigraphique. Le 
module volumique mailie comprend ainsi iui nombre de plan superieur au nombre de 
surfaces critiques. La mise en oeuvre d'un tel modele volumique implique aussi un 
10 travail consequent. Un modfele volumique peut presenter de 200 000 a un million de 
points, sur une trentaine de plans successifs. La fourniture d'une grille peut 
repr6senter un travail de 1 & 2 mois/homme. 

U existe done un besoin d'une solution qui permette de fournir simplement, 
pour des realisations obtenues a partir des mesures geophysiques, des grilles, 
15 autrement dit diss modules maill6s des realisations. 

Dans un mode de realisation, I'invention propose en consequence un procede 
de calciil d'une description mailiee d'une realisation d'un reservoir, la realisation 
comprenant unc pluralite des surfaces stratigraphiques, le precede comprenant : 

- la fourniture 

20 - d'une realisation de reference du reservoir, la realisation de reference 

comprenant une pluralite de surfaces stratigraphiques; 

d'une description mailiee de reference pour la realisation de reference; la 
description de reference comprenant une pluralite de plans dont certains 
d6crivent les surfaces stratigraphiques chaque plan comportant une pluralite 
25 de points; et 

d'au moins deux surfaces stratigraphiques de la realisation correspondant a 
deux surfaces stratigraphiques de la realisation de reference; 

- pour deux points homologues des deux plans d6crivant les deux surfaces 
stratigraphiques de la realisation de reference 

30 - le calcul des points sous-jacents des deux surfaces stratigraphiques de la 

realisation de reference; 
- le calcul des deplacements des points sous-jacents dans le passage des deux 

surfaces stratigraphiques de la realisation de reference aux deux surfaces 

stratigraphiques correspondantes de la realisation; les deux points sous- 
35 jacents deplaces constituant deux points homologues des plans de la 

description mailiee d6crivant les deux surfaces stratigraphiques de la 

realisation; 
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la foumiture de plans de la description maillee par interpolation entre les 
points homologues de ces deux plans. 
Dans un mode de mise en oeuvre, les deux surfaces stratigraphiques de la 
realisation comprennent la surface de toit et/ou la surface de base de la realisation. 
5 L'interpolation peut aussi comprendre, pour une autre surface stratigraphique 

de la realisation et pour deux points homologues des plans de la description maillee;, 
le calcul de 1'intersection entre l'autre surface stratigraphique et une droite passant par 
les deux points homologues, Intersection constituant un point du plan de la 
description maill6e decrivant l'autre surface stratigraphique. 
10 ^interpolation peut aussi comprendre, pour deux points homologues de deux 

plans de la description mailtee decrivant des surfaces stratigraphiques adjacentes, la 
foumiture d'un point d'un autre plan sur la droite reliant ces deux points homologues, 
dans les memes proportions de longueur que dans la description maillee de reference. 
Dans le cas oil les plans de la description maillee de reference s'etendent entre 
15 deux plans decrivant des surfaces stratigraphiques adjacentes sont paralleles, 
Interpolation peut comprendre, pour deux points homologues de ces deux plans, la 
foumiture d'un point d'un autre plan sur la droite reliant ces deux points homologues 
de sorte que les plans de la description maillee correspondant aux plans paralleles de 
la description maillee de reference soient paralleles. 
20 La description maillee de reference peut aussi comprendre au moins un point 

irr6gulier avec au moins deux positions; dans ce cas, le calcul d'un point irregulier 
correspondant de la description maillee comprend : 

la recherche de points reguliers de la description maillee de reference voisins 
d'une posit) on du point irregulier ; 
25 - le calcul d'une position du point irregulier correspondant dans la description 
maillee a piutir des deplacements des points r6guliers voisins entre la description 
maillee de reference et la description maillee. 

L'invention propose encore un programme de calcul d'une description maillee 
d'une realisation d'un reservoir, la realisation comprenant une pluralite de surfaces 
30 stratigraphiques, le programme comprenant : 
- une routine d 1 introduction 

d'une realisation de reference du reservoir, la realisation de reference 
comprenant une pluralite de surfaces stratigraphiques; 

d'une description maillee de reference pour la realisation de r6f6rence; la 
35 description de reference comprenant une pluralite de plans dont certains 

decrivent les surfaces stratigraphiques, chaque plan comportant une 
pluralite de points, et 
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d'au moins deux surfaces strati graphiques de la realisation correspondant a 
deux surfaces stratigraphiques de la realisation de reference; 

- une routine de calcul, pour deux points homologues des deux plans decrivant les 
deux surfaces stratigraphiques de la realisation de reference 

5 - des points sous-jacents des deux surfaces stratigraphiques de la realisation 

de reference; 

des displacements des points sous-jacents dans le passage des deux surfaces 
stratijjraphiques de la realisation de reference aux deux surfaces 
stratijjraphiques correspondantes de la realisation; 
10 les deux points sous-jacents deplaces constituant deux points homologues 

des plans de la description maillee decrivant les deux surfaces 
stratigraphiques de la realisation ; 

- une routine de calcul de plans de la description maillee par interpolation entre les 
points homolof naes de ces deux plans. 

15 Les routines peuvent anssi etre adapt6es a la mise en cxuvre des difffcrentes 

variantes du precede, 

D'autres earacteristiques et avantages de l'invention apparaitront k la lecture de 
la description detaillee qui suit des modes de realisation de l'invention, donn6s k titre 
d'exemple uniquement et en references aux dessins qui montrent : 
20 - figure 1, une representation schematique d'un reservoir avec les surfaces d'une 
realisation; 

figure 2 y une representation schematique en coupe dans un plan vertical des 
surfaces de deux realisations et des plans d'tuae grille de reference correspondant 
k Tune des realisations; 

25 - figures 3 a 6, un ordinogramme des etapes d'un procede selon un premier mode 
de mise en ^uvre de l'invention; 

- figure 7, une representation schematique analogue a celle de la figure 2, montrant 
la grille calouiee selon le procede des figures 3 a 6; 

figure 8, une representation schematique d'un nceud de faille dans une description 
30 maillee; 

figure 9 ? une representation des etapes d'un procede applique aux naends de faille 
dans un mode de realisation de l'invention. 
L'invention propose de determiner la grille ou modele maille d'une realisation a 
partir de la grille ou modele maill£ d&jk determine pour les surfaces d'une autre 
35 realisation. Elk repose sur la constatation que les differentes realisations derivent des 
mernes donn6e$ geophysiques et que les similitudes entre de$ realisations sont 
sufBsantes pour eviter de calculer k nouveau 1'ensemble de la grille pour une 
nouvelle realisation. 
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Les notations suivantes sont utilisees dans la suite de la description. On note 
S°i, S°2> ...» S°„ les surfaces d'une realisation R°, dite realisation de r<5f6reoce, pour 
laquelle on dispose d'une grille. Dans cette realisation de reference, S°i est la surface 
de toit, S° n est la surface de base et S°; ? 1 < i < n, sont les surfaces intermediaires. 
5 Comme indique plus haut, n a une valeur typiquement comprise entre 3 et 7 — la 
realisation presentant de 1 a 5 surfaces intermediaires. On appelle dans la suite 
"couche" le volume s'etendant entre deux surfaces de la realisation. Chaque surface 
S°i peut etre representee par un ensemble de points M°jj. 

On considere par exemple un rep6re dans lequel la direction z est la verticale 

10 orientee vers bas, les directions x et y fournissant Thorizontale. Le nombre de 
points pour une surface peut typiquement atteindre plusieurs mitliers. Comme 
explique plus haut, a chaque point peut etre associee une ou plusieurs incertitudes, 
habituellement exprim6es sous forme d'un vecteur. La representation des surfaces de 
la realisation de reference comme des autres realisations est sans incidence sur le 

15 precede decrit. On considfere £ titre d'exemple dans la suite une representation des 
surfaces S°i par des points M°ij. 

La grille G° associee a la realisation de reference est notee GV G\ G°m et 
est qualifiee dans la suite de grille de r6f<§rence. Dans cette notation G 0 , designe un 
plan de la grille, Un tel plan peut s'appuyer sur une surface S°i de la realisation; il 

20 peut aussi s'agir d'un plan s'<§tendant dans une couche entre deux surfaces de la 
realisation. L'utilisation du terme "plan" n'implique en aucune fa?on que les "plans" 
G°j pr6sentent ta structure g6om6trique d'un plan au sens classique de ce terme en 
mathematiques; en revanche, les "plans" G°i prdsentent la structure topologique du 
plan au sens classique. La grille de reference est formee d'un ensemble de points 

25 N° g .hj de coordonnees (u g , Vh, wO- Les points (u g , Vh> w$) pour une valeur fixe de i 
repr6sentant la mod61isation de la surface S°j. Les indices g et h sont des indices 
entiers dont les valeurs sont fonction du nombre de points dans la grille pour le plan 
G°i. On comprend de cette notation que la grille pr<6sente pour chaque plan une 
topologie similaire. On notera que les points de plans voisins peuvent Stre confondus, 

30 par exemple lorsque les surfaces sous-jacentes sont en biseau et qu'une couche 
disparait. La similitude entre les topologies des plans n'exclut done pas une vari6t6 de 
la g6om6trie de ces plans. 

Cette hypothese d f homologie dans la modelisation simplifie les caJculs 7 mais 
n'est en rien indispensable a Ja mi.se en oeuvre de Tinvention. Deux points (ug, Vh, Wji) 

35 et (u gJ Vh, Wi^), pour deux valeurs il et i2 differentes peuvent Stre qualifiees 
d'homologues dans les surfaces S°n et S°i2; autrement dit, du point de vue 
g6ologique, ils correspondent a une sedimentation sur un raeme point. On appelle 
alors "pilier" un ensemble de points (u^ v h> w,), pour des valeurs g et h fixes et une 
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valeur de i variant de 1 k m. Le pilier est note P gih . Du point de vue qualitatif, un 
pilier est constitue d'un ensemble de points superposes suivant la direction de 
sedimentation — avant deformation eventuelle. Cette hypothese sur la grille n'est pas 
indispensable a la realisation de l'invention, comme decrit dans la suite; il s'agit 
5 toutefois d'un mode de realisation correspondant h une grille foumie selon les 
solutions de J'etat de la technique. 

On note S k u S k 2 , S\ les surfaces d'une autre realisation R k . La realisation 
R k est representee par des points, tout comme la realisation de reference; dans cette 
premiere representation, la surface S k est representee par un ensemble de points 

10 M k j, correspondant aux points M°jj. On peut aussi representer la realisation par 
Tensemble des points M°jj de la realisation de reference, assorties d'un tirage des 
deplacements de chacun de ces points pour atteindre Vautre realisation; cette 
deuxieme representation de Tautre realisation correspond en fait k des valeurs 
correi6es des incertitudes aux diff6rents points. Dans une representation comme dans 

15 Tautre, il est possible de definir une homologie entre les points des deux realisations : 
deux points des deux realisations sont homologues dans la mesure ou ils representent 
le mgme point de la structure geologique sous-jacente, dans deux hypotheses 
correspondant «iux deux realisations. Pour la description de Tinvention, on considere 
dans la suite la premiere representation. 

20 Comme explique plus haut, le probleme de l'invention est de disposer d r une 

grille ou model e mailie pour cette realisation R k . L'invention propose de fournir cette 
grille a partir de la grille de reference C^. 

La figure 1 montre une representation schematique d'un reservoir avec les 
surfaces de la realisation de reference R°. On a repr6sent6 sur la figure les axes x, y et 

25 z du repere, ainsi que les surfaces S°i> S° 2 , S° n de la realisation R°, dans un 
exemple avec le nombre n de surfaces iimite a 3 pour des raisons de clarte de la 
figure. On a en outre repr6sent6 sur la figure 1 le maillage G°i du plan s'appuyant sur 
la surface de toit S°i de la realisation; du fait que le plan G°i s'appuie sur la surface 
S°i, celle-ci est difficilement visible sur la figure. La figure montre encore un plan de 

30 faille F traversant le reservoir; des details sur les failles sont donn6s plus bas. Enfin, 
on a represente. mais uniquement sur le bord du reservoir situe sur Tavant dans la vue 
en perspective de la figure I, les plans G°i & G 0 ^, en tra?ant les piliers 
correspondants. Pour la clarte de la figure, on n'a pas represente les autres points de 
ces plans. 

35 La figure 2 montre une representation schematique en coupe dans un plan 

vertical des sur faces de deux realisations et des plans d'une grille de reference; dans 
une situation roelle, il n'y a aucune raison pour que les points des surfaces ou les 
points de la grille soient dans un plan comme celui de la figure 2. La figure 2 n'est 

R:\Brrvets\2G200\20507EP.doi - 56/0S/U2 . 17:03 - 1/2* 



XVI • U • 



3 * JLU 



8 

done pas representative de ce point de vue d'une situation reelle, mats n'a pour but 
que de permettre Pexplication- 

La figure 2 montre un exemple simple d'une realisation de reference R° 
presentant deux couches definies par les surfaces S°t, S°2 et SV Ces surfaces sont 
5 representees en gras sur la figure 2. La grille de reference G° presente neuf plans G°i 
£ G° 9 . Les plant; G°i, G% et G° 9 "s'appuient" sut les surfaces S°i, S° 2 et S°3, autrement 
dit sont des modSlisations de ces surfaces; ils ne sont pas montr£$ sur la figure pour 
des raisons de clarte, dans la mesure ou ils sont tres proches des surfaces sur 
lesquelles ils s'appuient. Les plans G°2 & G°5 sont des plans intermediaires, dans la 

10 couche definie entre les surfaces S°i et S U 2; dans cette couche, les plans sont 
paralieies k la surface inferieure de la couche, ce qui correspond a une hypothese de 
sedimentation sur la couche S° 2 et d'erosion par la couche S°i. Les plans G° 7 et G° 8 
sont des plans intermediaires, dans la couche definie entre les surfaces S°2 et S 0 ^; 
dans cette couche, les plans sont repartis de maniere proportionnelle; ce qui 

15 correspond & une subsidance pendant la sedimentation. Autrement dit, les points d'un 
pilier sont rdgutierement repartis dans la couche. La figure montre aussi la realisation 
R k ; elle presente aussi trois surfaces S k i, S k 2 et S k 3, representees en pointilles sur la 
figure 2. Ces surfaces sont proches des surfaces correspondantes de la realisation de 
reference. 

20 La figure montre encore un pilier P° &h de la grille de reference. Ce pilier est 

une droite, qui dans Texemple passe dans le plan de la figure. On n T a reference sur la 
figure que les points N° g ,h,2 et N° gf h,8, pour plus de clarte. 

Comme cxplique plus haut, le probleme de Tinvention est de calculer ou definir 
k partir de la realisation de reference R°, de la grille de reference G° une grille pour la 
25 realisation G k . Le precede qui permet de definir cette grille est expliqu6 en reference 
aux figures 3 k 6, en utilisant les elements de la figure 2. 

Les figures 3 a 6 montrent un ordinogramme des etapes d'un proc6d6 selon un 
premier mode de mise en ceuvre de Tinvention, Dans ce mode de realisation, les 
piliers P° g ,h sont des droites. Cette caracteristique correspond a une contrainte de 
30 construction de la grille, qui est utile pour permettre Vutilisation de logiciels de 
* calculs d'ecouk ment de fluide dans le reservoir. 

Sur la figure 3, a la premiere etape 2, on fournit la realisation de reference R°, 
la grille de reference G° associ6e k cette realisation, ainsi que la realisation R k pour 
laquelle on souhaite obtenir une nouvelle grille. Cette etape rappelle les elements de 
35 depart du procede. 

A retape suivante 4, on considere un pilier P 0 ^ de la grille R°; par exemple on 
peut considerer lors de la premiere iteration le pilier P°u. 
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Aux etap^s 6 et 8, on determine un deplacement de ce pilier, lorsque Ton passe 
de la realisation R° k la realisation R\ On peut pour cela considerer le deplacement 
des intersections du pilier avec les surfaces de toit S°i et de base S° n lors du passage 
de la realisation R° a la realisation R k . Les intersections du pilier avec les surfaces de 
5 toit S°i et de base S° n sont en fait les points de ces surfaces sous-jacents aux 
extremites du pilier, autrement dit les points des surfaces qui sont represents par les 
extremites du pilier dans la realisation R°. On pourrait utiliser une autre methode 
pour determiner les points sous-jacents, par exemple en consid6rant le point de la 
surface le plus proche d'une extremite du pilier. 

10 A l'etape 6, on determine Intersection du pilier avec les surfaces de toit S°i et 

de base S°t, dims la representation de reference. On peut pour cela simplement 
considerer les extremites du pilier, c*est-a-dire les points N° gy h,i et N° s ,h,n et calculer 
les intersections I° B4 h,i et I° s ,h,n de la droite passant par ces deux points avec les 
surfaces S°i et S° n . Si Ton utilise pour la representation de la realisation R° un 

15 maillage avec des faces triangulares, ceci revient a determiner Intersection de la 
droite passant par les points N° &h a et N 0 ^™ avec une des faces de la surface. Le 
point l° gf )u s'exprime comme un barycentre des trois sommets M°i r i» M°i jT ^ et M 0 i, r 3 
de la face triangulaire de la surface S°i travers6e par le pilier, avec des coefficients 
oi\ f Pi et yi. II en est de meme pour le point l\h.ir», avec des coefficients cc^, p m et y m . 

20 Bien entendu, si les extremites du pilier sont aussi des points de la surface, il suffit de 
considerer les points en cause - qui reste un barycentre avec des coefficients 
particuliers. Comme indique plus baut, on n'a pas, pour plus de clarte, porte a la 
figure 2 les points l\ h ,i et l\h,9, ni les points correspondants de la surface. 

A retape 8, on determine le deplacement des points l°g.h,i et I 0 ^ lors du 

25 passage de la realisation R° a la realisation R k . En d'autres termes, on considere les 
points de la realisation R k qui correspondent aux points I°g,h,t et I^h.n ou sont 
homologues de ces points. Si les surfaces S k de la realisation R k sont representees 
par des points M\ r , comme indique plus haut, le point de la surface de toit peut 
simplement etre obtenu en considerant le barycentre des points M\h, M k u2 et M k i,r3 

30 affecte des mSmes coefficients aj 7 pi et y x . On precede de la meme fa9on pour la 
surface de base, avec les coefficients correspondants. Dans cet exemple, on utilise les 
coordonn6e$ barycentriques par rapport aux points les plus proches de la surface. 

A Tissue de cette etape, on trouve done des points l\ M et I k gi h, n des surfaces de 
toit S k i et S k n de base de Tautre realisation R k . Ces points constituent les extr6mit6s 

35 du pilier F\ h dans la grille R k . 

A retape 10, on passe au pilier suivant, et on recommence le precede pour ce 
pilier suivant, en bouclant vers retape 6. La figure ne montre pas explicitement la 

X:\BreveM\2020O\202O7EP. <Jv,- - 2$/0«/O2 - 17:08 . 9'24 



10 

sortie de la boucle, lorsque tous les piliers ont ete trait6s; elle ne montre pas non plus 
la fa9on de paroourir l'ensemble des piliers, qui peut etre quelconque. 

Apres avoir parcouru la boucle pour l'ensemble des piliers, on a defini les 
points N k g )M et N k 8>h>m formant les plans G k , et G k m s'appuyant sur les surfaces de toit 
5 et de base de 1'autre realisation. On a done determine simplement k partir 

- des surface;; S°i et S° n de toit et de base de la realisation de reference, 
des plans G°\ et G° n s'appuyant sur ces surfaces et 

- des surface,*; S k i et S k n de toit et de base de 1'autre realisation, 

les plans G k i zt G k m s'appuyant sur les surfaces de toit et de la base de 1'autre 
10 realisation R k . 

Aux etapes 12 k 58, on determine k partir des plans de la grille s T appuyant sur 
les surfaces de toit et de base de 1'autre realisation R k les plans intermediates de la 
grille R\ La figure 3 et 4 ainsi que les figures 5 et 6 montrent differentes solutions 
possibles pour interpoler les plans intermediaires. La premiere solution, representee 

15 sur la partie inferieuxe de la figure 3, s'applique aux plans intermediaires s'appuyant 
sur des surfaces de la realisation R k . La deuxifeme solution, representee sur la figure 
5, s T applique aux plans definis a l'interieur d'une couche, en mode proportionnel, 
comme dans la couche entre les surfaces S°2 et S°3 sur la figure 2. Ceci coiTespond a 
une tectonique synsedimentaire; la troisieme solution, representee sur la figure 6, 

20 s'applique aux plans definis a Tinterieur d'une couche, lorsque Ton souhaite que la 
grille G k reprocluise une hypothese g^ologique particuli^re presente dans la grille G° : 
il peut s'agir comme sur la couche entre les surfaces S°t et S°2 sur la figure 2 d'une 
Erosion (mode "On Lap 1 ' : erosion) ou encore de plans paralldles a la surface de toit, 
dans une hypothese dite "Top Lap", e'est-a-dire transgression-regression. 

25 La partie inferieure de la figure 3, etapes 12 k 22, permet de definir les plans de 

la grille G k s'appuyant sur les surfaces intermediaires de la realisation. A 1'etape 12, 
on considere une surface intermediaire S k de 1'autre realisation. Cette surface est 
hoirtologue d'une surface S°i de la realisation de reference, sur laquelle s'appuie un 
plan G°j de la grille de reference. On considere k l'6tape 14 un pilier donne, par 

30 exemple Ie pilier P k g . h de la grille R k : ce pilier est defini comme la droite - ou le 
segment de droite si les surfaces consid6r6es aux etapes 4 a 10 sont la base et le toit - 
passant par deux points homologues N k g ,h,i et N k gi h,„. 

A retape 16, on calcule Intersection de ce pilier avec la surface intermediaire 
S k . On peut pioceder comme explique plus haut a 1'etape 6. Cette intersection est 

35 consid6ree corr me le point N k g hj du plan G k s'appuyant sur la surface intermediaire 

s\ 

A Tetape 18, on passe au pilier suivant, et on recommence le calcul pour ce 
pilier suivant, en bouclant vers Tetape 16. S'il n'y a plus de pilier suivant, on passe k 
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T6tape 20. A l'etape 20, on considere la surface suivante et on boucle vers T6tape 14. 
S'il n'y a plus d<3 surface suivante, on passe a l f 6tape 22. 

Ces 6tapcs fournissent sirtiplement pour la surface interm6diaire S k Ie plan 
correspondant G k correspondant. Ce precede preserve la propriete qu'un pilier de la 
5 description maillee G k est une droite. Cette propriete peut etre utile, notamment pour 
permettre d'utiliser ensuite des logiciels d^coulement de flux dans le reservoir. Elle 
preserve aussi une homologie entre les plans s*appuyant sur les diff6rentes surfaces, 
autreinent dit uae similitude entre les topologies des surfaces. Le precede des etapes 
12 k 22 pr6serv 2 aussi la rectitude des piliers. 

10 II est manifeste que le precede des etapes 12 a 22 ne fonctionne que dans la 

mesure oil l'6tape 16 est possible - c'est-a-dire dans la mesure oii il est possible de 
calculer Pinters ection de la surface S k avec le pilier. Le precede des etapes 12 k 22 
est done adapts pour foumir les plans de la grille G k s'appuyant sur les surfaces de 
cette realisation . Dans la mesure ou ces plans sont obtenus par intersection des piliers 

15 avec les surfaces, ils fournissent un module mail!6 pr6cis des surfaces* 

II reste done a calculer ou definir les plans intermediaires de la grille G k , qui ne 
s'appuient pas sur des surfaces interni6diaires de la realisation R k . Le precede 
propose a cet effet deux solutions; il est avantageux, comme expliqu£ ci-dessous, de 
mettre en oeuvoi des deux solutions, en fonction des couches rencontrees. On pourrait 

20 aussi ne mettre en oeuvre qu'une seule de ces deux solutions pour Tensemble de la 
realisation. 

Sur la figure 4, k l'6iape 24, on considfere une couche — autrement dit deux 
surfaces adjacentes de la realisation R k , pom: lesquelles on dispose d6ja des plans 
correspondants de la grille, d6flnis aux Stapes 12 k 22. On peut consid6rer lors de la 
25 premiere iteration la premiere couche, entre les surfaces S°i et S°2. 

A l'£tape 26 suivante, on determine le mode d'interpolation des plans de la 
couche. Dans le cas d'une interpolation gSologique, on passe a l'etape 30 et dans le 
cas d'une interpolation homoth6tique, on passe k l'etape 28. La definition respective 
des plans de la couche aux etapes 28 et 30 est representee sur les figures 5 et 6. 
30 A l'6tape 32, les plans de la couche ont 6t6 d6finis. S'il existe encore une 

couche a traiter, on passe a la couche suivante, en bouclant vers l'etape 26. Sinon, le 
proc6d6 est termin6. 

Comme indique plus haut, Tinteret de rordinogramme de la figure 4 est de 
proposer des solutions alternatives pour le calcul par interpolation des plans 
35 intermediaires. On pourrait aussi ne choisir qu'une solution appliquee a toutes les 
couches definies entre les surfaces de la realisation R k . 

La figure 5 montre une solution d'interpolation des plans intermediaires, dans 
une hypothese ou les plans intermediaires de la realisation R k dans la couche en 
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cause sont homothetiques des plans interm6di aires de Ja couche en cause dans la 
realisation dc r£f6rence R°. La figure 5 correspond & T6tape 28 de la figure 4. La 
solution de la figure 5 repose sur la definition des plans G k de la description maill6e 
G k dans les mSmes proportions que dans la description mailtee de reference G°. 
5 Autrement dit, on importe dans la description mail!6e G k les positions relatives des 
plans G°i de ki couche en cause dans la description maillee de reference G°. La 
solution des etapes 36 a 42 est decrite pour une couche, entre deux plans G k p et G k q 
s'appuyant sur deux surfaces S k et S k +1 de la realisation R k . Dans l'exemple de la 
figure 2, on considere la couche s'etendant entre les surfaces S k 2 et S k 3, sur lesquelles 

1 0 s'appuient les p lans G k 6 et G%. 

On considere a l'etape 36 un pilier donne, par exemple le pilier P 1 ^ de la grille 
R k : ce pilier est defini dans la couche consid6ree comme le segment de droite 
s'etendant entre les points N k g ,h. p et N k gi h,q des plans G k p et G k q . A la premiere 
iteration de la boucle de la figure 4, on peut simplement considSrer le pilier P k i,i- 

15 A l'etape 38, on calcule, pour le pilier correspondant P° g ,h dans la description 

maill£e de r6f6rence G°, la position des points N°g,h,t, p < t < q, des plans de la grille 
de reference s'etendant a Tinterieur de la couche en cause dans la grille de reference. 
La position des diff&rents points est donn6e par exemple par le rapport des longueurs 
aux extremites du segment, ou, ce qui est equivalent, par les coefficients du 

20 barycentre par rapport a ces extremites. On pourrait aussi utiliser des coordonnees 
dans un rep^re orthonorm6 dont le vecteur unitaire relie les extr6mit6s du pilier* 

A Tetape 40, on definit les points N k g ,h.t sur le pilier P k &h, de sorte a respecter 
sur ce pilier les mSmes proportions que sur le pilier P° g , h . Ceci s'effectue par un 
calcul des rapports de longueurs sur le segment [N k g,h, p5 N k gt h. q ], °u encore par un 

25 calcul du barycentre des extr6mitds du segment ou enfin par le calcul d'un vecteur. 
Cette etape pennet de definir les points des differents plans de la couche sur le pilier. 
Dans l'exemple de la figure 2, on ddfinit de la sorte les points N k g ^ 4 7 et N k g ,h.s> ces 
points divisant le segment [N k g .h,6» N k gj h,9] dans les m6me$ proportions que les points 
N° g> h,7 et N° g ,h t 8 divisent le segment [N° s ,h,6, N\h,<>] dans la grille de rdfdrence. 

30 On passe ensuite au pilier suivant, a Tetape 42, en bouclant vers l'etape 38. S'il 

n'existe pas de pilier suivant, on retourne & T6tape 32 de la figure 3. 

On arrive done a definir les plans de la description maillee G k dans une couche 
entre les plans s'appuyant sur deux surfaces successives. Dans cette couche, les 
positions relatives des plans sont les memes que celles des plans correspondants dans 

35 la grille de reference. Ainsi, dans l'exemple de la figure 2, les plans G° 7 et G°s 
divisent chaque pilier en trois parties sensiblement egales dans la couche s'etendant 
dans la grille de reference entre les surfaces S°2 et S°3; il en est de meme des plans 
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G k 7 et G k g dans la couche s'etendant dans la grille G k entre les surfaces S k 2 et S k 3 ou 
les plans G k 6 e1 GV Ceci est repr6sent6 k la figure 7. 

Comme dans le cas de la figure 3, Interpolation proposee k la figure 6 
preserve lTiomologie des topologies des plans intermedi aires; elle preserve aussi le 
5 caractere rectili gne des piliers. 

La figure 6 montre une solution ^interpolation des plans intermediaires, dans 
une hypothese ou les plans intermediaires de la r6a!i$ation R k dans la couche en 
cause reproduisent la meme hypothese geologique que les plans intermediaires de la 
couche en cause dans la realisation de reference R°. Cette figure correspond k retape 
10 30 de la figure 4. On considere dans l'exemple de la figure 6 deux hypotheses 
geologiques, a savoir 

des plans ititerm6diaire$ paralleles au plan inferieur de la couche en cause; ou 
des plans ir termediaires paralleles au plan sup6rieur de la couche en cause. 
La premiere h;/poth£$e correspond a une sedimentation sur la surface sur laquelle 
15 s'appuie le plan inferieur et le cas 6ch6ant a une erosion en surface. La deuxieme 
hypothese correspond a une sedimentation avec une regression ou transgression 
marine. Sur la figure 2, la couche entre les surfaces S°i et S°2 presente des plans 
paralleles G°2 ii G°s paralleles au plan inferieur G°6 s'appuyant sur la surface S° 2 ; la 
repartition des plans dans la couche correspond a la premiere hypothese geologique 
20 proposee plus Itaut. 

La solution de la figure 6 repose done sur la definition des plans G k de la 
description maillee G K de 1'autre realisation R k , de sorte a respecter aussi dans la 
couche en cause les memes hypotheses geologiques que dans la description maillee 
de reference G°. Autrement dit, pour les hypotheses geologiques de paralieiisme 
25 considers, on importe dans la description maillee G k les contraintes de paralieiisme 
des plans G°\ d« la couche en cause dans la description maillee de reference G°. 

La solution des etapes 36 a 42 est decrite pour une couche, entre deux plans 
G k p et G k q s'aspuyant sur deux surfaces S k et S k i+i de la realisation R k . Dans 
l'exemple de la figure 2 y on considere la couche s'etendant entre les surfaces S k \ et 
30 S k 2, sur lesqueDes s'appuient les plans G k i et G k $. 

On considere k Vetape 50 une hypothese geologique; dans l'exemple considere 
ceci revient a ohoisir celui des plans sup6rieur et inferieur de la couche auquel les 
plans intermediaires sont paralleles. 

A t'etape suivante 52, on determine le pilier le plus long de la couche. Ceci 
35 revient k maximiser sur Tensemble des piliers la longueur du segment [N k g ,h.p, 
N k g ,h.qL On trouve le pilier P k s , h le plus long. 

On interpole ensuite les plans intermediaires sur le pilier le plus long. On peut 
simplement repartir sur le pilier le plus long de fa?on regulifere les points N k g .h.i ? p < t 
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< q, des plans intermediates; alternativement, comaie represent^ a Tetape 54, on 
peut rechercher la position des points N° g ,ii,t, p < t < q sur le pilier de la grille de 
reference P° gv h oorrespondant au pilier le plus long. A T6tape 56, on d6finit les points 
N k & h,t sur le pilier P k g .h, de sorte a respecter sur ce pilier les memes proportions que 
5 sur le pilier P° g ,h Ceci peut s'effectuer comme explique en reference a Tetape 40 de la 
figure 5. 

A l'etape 58, on calcule les points N k g \h\i sur les autres piliers P k s ^» avec (g\ h') 
different de (g : h), de sorte & respecter Thypoth^se de parallelisme. II suffit, dans 
Texemple d ! un parallelisme avec le plan s'appuyant sur la surface inferieure de la 
10 couche, que les distances entre un point d'un plan intermediate et le point du mSme 
pilier du plan inf&rieur soient 6gales pour tons les piliers, soit 

d(N k g ..h> M k g^.q) = d(N k gfh .c N k g , h>q ) 
pour tout triplet (t, g', h ? ) v6rifiant p < t < q et (g', h') different de (g, h). II reste 
possible a cette etape qu'un plan soit confondu partiellement avec le plan de base ou 
15 de toit — par exemple si la longueur d'un pilier devient trop faible par rapport a la 
longueur du pilier le plus long. Le precede de la figure 6 permet done aussi de 
repr6$enter des plans terminant en biseau. 

Apres l'etape 58, on retoume a l'etape 32 de la figure 3. 

On arrive done 4 d6fimr les plans de la description maillee G k dans une couche 
20 entre les plans s'appuyant sur deux surfaces successives. Dans cette couche, le 
parallelisme des plans est identique au parallelisme des plans correspondants dans la 
grille de r6f6rence. Ainsi, dans Texemple de la figure 2, les plans G°2 i G° 5 sont 
paralleles au plan G°6 qui s'appuie sur la surface S°2. De meme, les plans G k 2 a G k 5 
sont paralleles au plan G°$ dans la grille G k . 
25 La figure 7 montre une representation analogue a celle de la figure 2, sur 

laquelle on a repr6sent6 

en traits po;ntilles, les trois surfaces S°i, S°2 et S°3 de la realisation de reference 
R°, ainsi que le pilier P° g ,h; 

les trois surfaces S k i, S k 2 et S k 3 de la realisation R k ainsi que la grille G k obtenue 
30 par le proc6d6 d6crit en reference aux figures 3 4 6, avec le pilier P k 5 ,h. 

On constate sur la figure que les plans intermediates dans la couche superieure 
sont comme sur la figure 2 paralldles 4 la surface inferieure de la couche, tandis que 
les plans intermediaires dans la couche inferieure sont regulierement repartis dans 
cette couche. Encore une fois, la figure 2 comme la figure 7 correspondent a des 
35 exemples simples, dans lesquels les surfaces sont sensiblement planes et 
sensiblement piiralteles: le proc6d6 s'applique aussi pour des surfaces gauches ou des 
formes complexes. 
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Le procede decrit aux figures 3 a 6 perniet ainsi de calculer une description 
maillee G k d'une realisation R k , k partir 
- d'une realisation de reference R°; 

d'une description maillee de reference G° de cette realisation de reference; et 
5 - d'au moins deux surfaces stratigraphiques S k \ , S k n de la realisation R k . 

Le procede n'implique pas de calculer k nouveau toute la description maillee pour la 
nouvelle reali?5ation, mais suppose simplement la connaissance des points 
homologues des surfaces a rint6rieur d'une realisation. H suppose aussi un calcul de 
deplacement, c'est-a-dire une homologie de la representation entre les di verses 

10 realisations. Ceci est simplement assure dans les exemples proposes plus haut, si les 
realisations correspondent k des tirages d'incertitudes, ou encore si les realisations 
presentent une meme topologie. 

Le calcul d'une nouvelle description maillee G k suivant le precede de 
Tinvention peut s'effectuer en quelques minutes de temps de calcul. 

15 Le proc6de d6crit en reference aux figures 3 k 6 s'applique k toutes les surfaces 

et permet de fournir une description maillee. On decrit maintenant en reference aux 
figures 8 et 9 un mode de realisation dans lequel on applique un traitement distinct a 
des points des surfaces definissant des failles. Le mode de realisation de la figure 9 
n'est pas obligatoire, mais permet d'obtenir une description maillee plus continue ou 

20 moins perturbee par les failles. 

Le mode de realisation des figures 8 et 9 repose sur Texistence, dans certaines 
descriptions niaill6e$ de surfaces, de points particuliers. Ces points sont 
repr6sentatifs cle failles ou de discontinuites dans la surface decrite. Pour rendre 
compte de ces discontinuites, la description maillee peut utiliser des points 

25 particuliers, qualifies de "split nodes" ou noeuds de coupure. Un tel point de la 
mod61isation peut presenter plusieurs positions geomdtriques correspondant k deux 
ou plusieurs panneaux de faille. La figure 8 montre une representation schematique 
d'un nceud de fsiille dans one telle description maall6e, dans un exemple ou les points 
irreguliers presentent deux positions. On a represents sur la figure une surface et trois 

30 plans de part et d' autre de la surface. La surface sous-jacente presente une faille 
separant deux parties 62 et 64 de la surface; dans l'exemple, la faille est sensiblement 
verticale. Les points des plans autour de la partie 62 ou de la partie 64 sont dits 
r£guliers; ils sont representes dans 1'exemple de la figure 8 par des points de forme 
circulaire sur U>s plans. Un point de la description maillee au niveau de la faille est 

35 appele nceud ou point de coupure ou encore point irregulier; il presente deux 
positions geoxri6triques representees par des croix sur la figure 8; Tune de ces 
positions est utilisee pour definir le bord d'une partie de la surface, tandis que Tautre 
position est utilis£e pour definir l'autre partie de la surface. On a reference sur la 
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figure la premiere position 66 d'un point irregulier, qui est utilisee pour definir le 
bord d'un plan sur un cdt£ de la faille. L r autre position 68 du point irregulier est 
utilisee pour dSfinir le bord du plan sur l'autre c6t6 de la faille. Les deux positions du 
point irrdgulier sont sur la meme surface reglee 70 definissant le panneau de faille 
5 Le precede de la figure 9 propose de prendre en compte de fa?on particuliere 

les points irr6guliers de la surface. Pour chacune des positions g6om6triques des 
points irr^guliers, le deplacement est obtenu a partir du deplacement des points 
reguliers voisir*s. Cette solution permet, dans le modfele mailte obtenu, de mieux 
rendre compte des failles ou des irregularites. En particulier, cette solution assure que 

10 le panneau de faille reste represents par une surface reglee. 

Le procede de la figure 9 commence par un calcul des points r6guliers de la 
description maillee G k de la realisation R k , suivant le procede decrit en reference aux 
figures 3 i6. L'6tape 72 representee sur la figure 9 est done un calcul des points 
reguliers des descriptions raaillees des plans G k o et G k ra s'appuyant sur les surfaces de 

15 base et de toit S k o et $\ 9 puis un calcul par interpolation des points r6guljers des 
autres plans G k . 

A retape 74, on consid&re un point irregulier. Ce point est caract6ris£ par au 
moins deux positions, comme explique en reference a la figure 8. Soit done N° B ,h,t un 
point irregulier, qui dans Texemple fait partie du plan G°i s'appuyant sur la surface de 

20 toit et T 0 ^,! et U Q 5 ,h,i les deux positions de ce point. 

A Tetape 76, on determine, pour une position de ce point — par exemple pour la 
position T° Sth ,i lors de la premiere iteration - les points r6guliers voisins dans le plan 
G°i. On peut considerer les points reguliers les plus proches, ou considerer en outre 
d'autres points reguliers. On pourrait aussi considerer les positions voisines des 

25 points irregulkrs, pour lesquels le deplacement a 6te calcule lors des iterations 
precedentes. Le nombre de points considere a cette etape n'est pas essentiel - et on 
pourrait a la limite ne considerer qu'un point voisin. 

A retape 78, on calcule le displacement de la position du point irregulier — 
autrement dit la position T k g ,h,i dans le plan G k i 4 partir des d£placement$ des points 

30 reguliers voisins. On peut simplement effectuer une moyenne de ces deplacements, 
dans le cas le plus simple; on peut aussi ponderer les deplacements, par exemple en 
fonction des distances entre la position et les points reguliers voisins, ou suivant 
toute autre methode. Dans le cas le plus simple, on proedde k un simple calcul de la 
moyenne des deplacement et on trouve : 

-rO -rk _ Pip: * 

Plp2 



R:\Brcvcts\20200\20207EP.dra -26/03/02- 17:08- 16/24 



17 

avec N°pi #P 2,i les points reguliers voisins du point N° g ,h,i, le terme au denominates 
repr6sentant simplement le nombre de points r6guliers coixsideres. 
A Tissue de FStape 78, on obtient done la position d'un point irregulier dans la 
description mai ll£e GV 
5 A l'etape 80, on passe k la position suivante du point irregulier et on boucle 

vers retape 76 Si on a traite toutes les positions du point irregulier en cause, on 
passe a Tetape 32. 

A retape 82, on passe au point irregulier suivant et on boucle vers Tetape 74, Si 
on a traite tous les points irreguliers, le precede se termine. 

10 Les etapos 72 a 82 pennettent ainsi de g6n6rer les points irreguliers dans la 

description maillee G k de la realisation R k . Le traitement des positions des points 
irreguliers a partir des points reguliers voisins permet de preserver dans la 
description maillee G k les failles ou les autres irregularis exxstant dans la 
description maillee de reference. 

15 Le proc6d6 des figures 3 a 6, comme celui de la figure 9, permet de garder pour 

les plans G°k de la description maillee G k de toutes les realisations R k une topologie 
similaire A la topologie de la description maillee S° de la realisation de reference R p . 
On rend ainsi compte du fait que la description maillee est calcutee pour une 
realisation R k similaire a la realisation de reference : le procede prend en compte et 

20 exploite Thomologie qui existe entre les realisations. Dans la mesure oil la propriete 
de rectitude des piliers est assur6e pour les points reguliers, elle est aussi presence 
pour les points irreguliers dans le precede de la figure 9 — dans la mesure ou Ton 
considere, pour un pilier form6 de points irreguliers, les mernes choix des points 
voisins dans chaque plan, 

25 En particulier, le precede de la figure 9 assure une bonne representation des 

failles ou des irregularites; les failles representees par une surface r£giee dans la 
description maillee de reference G° sont aussi representees par une surface r6glee 
dans la description maillee G k obtenue par le procede. 

Les procides des figures 3 a 6 et 9 peuvent etre mis en ceuvre par des 

30 programmes appropries, avec des routines correspondent aux differentes etapes des 
precedes. La programmation de telles routines est A la portee de Thomme du metier, 
compte tenu des indications fournies plus haut en reference aux figures. 

Afin de rendre la programmation de Tinvention relativement aisee, il est 
preferable d'utiliser un langage de haut niveau, qui permette une programmation de 

35 type objet tel que le C++ ou le langage JAVA . Ce type de traitement est de 
preference inte^jre a un modeleur ayant une description a la fois surfacique et maillee 
afin de connaitre les relations topologiques entre les elements tels que points, piliers, 
couches et surfaces. En outre les logiciels de type modeleur decrivent deji les 
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fonctionnalit6s de base telles que les calculs d'intersection. Le logiciel vendu par la 
societe TSURF sous la marque GOCAD est un bon exemple de ce type 
d'environnement de programniation. 

Bien entendu, la pr6sente invention n'est pas limit6e aux modes de realisations 
5 d6crits k titre d'exemple. Dans les exemples proposes, on utilise pour 1'interpolation 
des piliers les surfaces de toit et de base de la realisation. On pourrait aussi utiliser 
deux autres surfaces de cette realisation. L'utilisation des surfaces de toit et de base 
presente Tavantage d'une distance maximale entre les surfaces, ce qui assure une 
precision optimale dans Interpolation du pilier, pour une incertitude donnee dans le 

10 plan de la surface. Si Ton utilise ces surfaces de toit et de base, on pent en toute 
rigueur parier de surfaces iutermedi aires, comine plus haut Meme si Ton n'utilise pas 
les surfaces de toit et de base, le proc6d6 s'applique aussi aux surfaces situees au- 
dessus et au-dessous des surfaces utilisees — bien qu'il ne s^gisse pas a proprement 
parler de surfaces "interm6di aires". 

15 Par ailleu rs, dans rordinogramme de la figure 3, le precede est explique avec 

un calcul pilier par pilier; on peut changer l'ordre des calculs proposes, par exemple 
en commen9ant par la recherche sur la surface de toit des points de passage des 
piliers, puis la recherche des points sur la surface de base et enfin interpolation des 
piliers. 

20 Uetape 6 de rordinogramme de la figure 3 suppose aussi que le pilier est une 

droite; dans la mesure ou le pilier n'est pas une droite, on peut consid6rer la 
projection des points N° g ,h,i et N 0 ^^ sur les surfaces S°i et S 0 n , de la meme mani&re 
que pour les calculs relatifs au toit et £ la base tels que ddcrit k P6tape 8 et non pas 
intersection proposee a Tetape 6 de la figure 3. Cette consideration s'applique 

25 naturellement pour les surfaces de faille qui ne sont pas n£cessairement des surfaces 
reglees. En tout etat de cause, on obtient les points et I° gj h, n de ]a droite passant 

par ces deux points avec les surfaces S°i et S° n . 

Comme pour Tordinogramme de la figure 3, on peut changer l'ordre des etapes 
de rordinogramme des autres figures, notamment pour ce qui est du balayage des 

30 surfaces et des points irreguliers. 

II n'est pas indispensable de mettre en oeuvre la solution de la figure 9 en 
Tabsence de points irreguliers. 

On rappelle encore que la moduli sation utilisee pour les surfaces dans les 
differentes realisations est sans incidence. Dans Texemple, on a considere des 

35 surfaces moduli s6es par des points; on a aussi consid6re une collection de points pour 
chaque realisation. On pourrait tout aussi bien mod61iser les differentes r6alisations 
en fournissant un ensemble de points pour la realisation de reference R° et n 
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ensembles de vecteurs d'incertitudes correspondents aux n autres realisations. On 
obtient de la m«kne fa£on une modelisation des surfaces des autres realisations. 

Le choix d'un mode d'interpolation dans une couche peut Stre une donnee de 
depart; il peut uussi etre deduit de l'analyse de la position des plans dans la grille de 
5 reference. 
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REVINDICATIONS 

1. Un procette de calcul d'une description mailiee (G k ) d'une realisation (R k ) d'un 
reservoir, la realisation comprenant une plurality des surfaces stratigraphiques (S k i, 

5 S k 2? . . S k n ), le precede comprenant : 

- la fourniture(2) 

d'une realisation de reference (R°) du reservoir, la realisation de reference 
comprenant une pluralite de surfaces stratigraphiques (S°i, S° 2 , . . •> S° n ); 

d'une description mailiee de reference (G°) pour la realisation de reference; la 
1 0 description de reference comprenant une pluralite de plans (G°i , G\ . >. > G° m ) dont 
certains d6crivent les surfaces stratigraphiques (S°i, S°2, . . S° n ), chaque plan (G°i) 
comportant une pluralite de points (N° g? h,i); et 

d'au moiris deux surfaces stratigraphiques (S k i, S k n ) de la realisation (R k ) 
correspondant ii deux surfaces stratigraphiques (S° 1? S° n ) de la realisation de 
1 5 reference; 

- pour deux poiuts homologues (N° g ,h,i» N° g ,hjn) des deux plans (G°i, G° m ) decrivant 
les deux surfaces stratigraphiques (S°i, S° n ) de la realisation de reference 

le calcul (6) des points (I° g ,h,i, I°g,h,Ti) sous-jacents des deux surfaces 
stratigraphiquer> (S°i, S° n ) de la realisation de reference; 

20 le calcul (8) des displacements des points sous-jacents dans le passage des deux 

surfaces stratigraphiques (S°i, S° n ) de la realisation de reference aux deux surfaces 
stratigraphiques correspondantes (S k i, S k n ) de la realisation (R fc ); 
les deux points sous-jacents deplaces constituant deux points homologues des plans 
(G k i, G k nO de Ui description mailiee (G k ) decrivant les deux surfaces stratigraphiques 

25 (S k i , S k n ) de la realisation (R^; 

- la fourniture de plans (G k ) de la description mailiee (G k ) par interpolation entre les 
points homologues de ces deux plans (G k i, G 1 *™)- 

2. Le precede de la revendication 1, caracteris6 en ce que les deux surfaces 
stratigraphiques (S°!, S° n ) de la realisation (R k ) comprennent la surface de toit et/ou 

30 la surface de base de la realisation. 

3. Le precede de la revendication 1 ou 2, caract6rise en ce que l'mterpolation 
comprend, pour une autre surface stratigraphique (S k j) de la realisation (R k ) et pour 
deux points homologues des plans (G k i, G k m ) de la description mailiee (G k ), le calcul 
(16) de l'intersection entre Tautre surface stratigraphique et une droite passant par les 

35 deux points homologues, Intersection constituant un point du plan (G k ) de la 
description mailiee (G k ) decrivant Tautre surface stratigraphique (S k ). 
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4. Le precede de la revendication 1 , 2 ou 3, caracterise en ce que Interpolation 
comprend, pour deux points homologues de deux plans (G k p> G k q ) de la description 
maillee (G k ) decrivant des surfaces strati graphiques (S k u S k +i) adjacentes, la 
fourniture (40) d'un point d ! un autre plan sur la droite reliant ces deux points 

5 homologues, dans les mSmes proportions de longueur que dans la description maillee 
de reference. 

5. Le proc6d6 de la revendication I, 2 ou 3, caract£ris£ en ce que des plans de la 
description maillee de reference s'etendant entre deux plans (G 0 P , G° q ) decrivant des 
surfaces strati graphiques adjacentes (S k > S k +t) sont paralleles, et en ce que 

10 Interpolation comprend, pour deux points homologues de ces deux plans (G k p , G k q ), 
la fourniture (54, 56) d'un point d r un autre plan sur la droite reliant ces deux points 
homologues de sorte que les plans de la description maillee correspondant aux plans 
paralleles de la description maillee de reference soient paralleles. 

6. Le proc6d6 de Tune des revendications 14 5, caracterise en ce que la 

15 description maillee de reference (G°) comprend au moins un point irregulier (N° E .h F i) 
avec au moins deux positions (T°g >hil> U°g,h.i) et en ce que le calcul d f un point 
irregulier correspondant de la description maillee (G^) comprend : 

- la recherche de points reguliers de la description maillee de reference voisins d'une 
position du point irregulier ; 

20 - le calcul d'une position du point irregulier correspondant dans la description 

maillee k partir des d6placements des points reguliers voisins entre la description 
maillee de reference (G°) et la description maillee (G k )- 

7. Un prognimme de calcul d'une description maillee (G k ) d'une realisation (R k ) 
d'un reservoir, ]a realisation comprenant une pluralite de surfaces stratigraphiques 

25 (S k i, S k 2 , . . S k n ), le programme comprenant : 

- une routine d^ntroduction 

d'une realisation de reference (R°) du reservoir, la realisation de reference 
comprenant une pluralite de surfaces stratigraphiques (S°i, S°2, S° n ); 

d'une description maillee de reference (G°) pour la realisation de reference; la 
30 description de i6f6rence comprenant une pluralitd de plans (G°u G°2, . . G° m ) dont 
certains decrivent les surfaces stratigraphiques (S°i, S°2, . S 0 *), chaque plan (G°0 
comportant une pluralite de points (N° 6 ,h,i); et 

d f au moins deux surfaces stratigraphiques (S k i> S k n ) de la realisation (R k ) 
correspondant a deux surfaces stratigraphiques (S°i, S° n ) de la realisation de 
35 reference; 
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- une routine de caicul, pour deux points homologues CN° s ,h > i, N\h.m) des deux plans 
(G°i 5 G 0 m ) decri vant les deux surfaces strati graphiques (S°i, S° n ) de la realisation de 
r^fSrence 

des points; (I° & .h,i ? I°g ? h,n) sous-jacents des deux surfaces stratigraphiques (S°j ? 
5 S° n ) de la realisation de reference; 

des deplaeements des points sous-jacents dans le passage des deux surfaces 
stratigraphiques. (S°t, S° n ) de la realisation de reference aux deux surfaces 
stratigraphiques. correspondantes (S k i, S k n ) de la realisation (R k ); 
les deux points sous-jacents deplac6s constituant deux points homologues des plans 
10 (G k u G k m) de la description maillee (G k ) decrivant les deux surfaces stratigraphiques 
(S k i, S k n ) de la realisation (R*); 

- une routine de caicul de plans (G k ) de la description maillee (G k ) par interpolation 
entre les points homologues de ces deux plans (G k i, G k ™). 

8. Le programme de la revendication 7, caracterise en ce que les deux surfaces 
15 stratigraphiques; de la realisation (R k ) comprennent la surface de toit et/ou la surface 

de base de la re alisation. 

9. Le programme de la revendication 7 ou 8, caract6ris6 en ce que la routine de 
caicul par interpolation est adaptee k effectuer, pour une autre surface stratigraphique 
(S k ) de la realisation (R k ) et pour deux points homologues des plans (G k i, G k m ) de la 

20 description maill6e (G k ), le caicul de Tintersection entre Tautre surface 

stratigraphique et une droite passant par les deux points homologues, Intersection 
constituant un point du plan (G k ) de la description maillee (G k ) ddcrivant 1'autre 
surface stratigraphique (S k ), 

1 0. Le programme de la revendication 7, 8 ou 9, caract6ris6 en ce que la routine de 
25 caicul par inter] >olation est adaptee a calculer, pour deux points homologues de deux 

plans (G k p , G k q ) de la description maillee (G k ) decrivant des surfaces stratigraphiques 
(S\ S k i+i) adjacentes, un point d'un autre plan sur la droite reliant ces deux points 
homologues, dans les meraes proportions de longueur que dans la description maillee 
de r6f6rence. 

30 11. Le programme de la revendication 7, 8 ou 9, caracterise en ce que des plans de 
la description maillee de reference s'etendant entre deux plans (G° p > G° q ) d6crivant 
des surfaces stratigraphiques adjacentes (S v u SVi) sont paralleles, et en ce que la 
routine de caicul par interpolation est adapt6e i calculer, pour deux points 
homologues de ces deux plans (G k p , G k q ), un point dW autre plan sur la droite reliant 
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ces deux points homologues de sorte que les plans de la description maill6e 
correspondant aux plans paralteles de la description maillee de reference soient 
paralteles. 

12, Le programme de Tune des revendi cations 7 a 11, caracterise en ce que la 
5 description maillSe de reference (G°) comprend au moins tui point irregulier (N° g .ha) 
avec au moins deux positions (T° & h,n U° fi .h.O et en ce que le programme comprend en 
outre une routine de calcul dun point irregulier correspondant de la description 
mailtee (G K ), la dite routine de calcul comprenant : 

- une sous-routine de recherche de points reguliers de la description maillee de 
10 reference voisins d'une position du point irregulier ; 

- le calcul d'un<5 position du point irregulier coirespondant dans la description 
maill6e k partir des deplacements des points reguliers voisins entre la description 
maillee de refeience (G°) et la description maillee (G k ), 
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ABREGE 

i 
i 

Uti proc6d6 de calcul d'une description mailtee (G k ) d\me realisation (R k ) d'un 
reservoir utilise j 
5 - une realisation de reference (R°) jdu reservoir, avec plusieurs surfaces 
stratigraphiques (S°i, S° 2 , . . S° n ); j 
- une description maill6e de ref6rence (<G°) pour la realisation de reference, avec 
plusieurs plans (G°i, G\ .... G°„) dont certains decrivent les surfaces 
stratigraphiques (S°i , S° 2 , . . . , S° n ); | 
10 - deux surfaces stratigraphiques (S k i, S k n ) de la realisation (R k ) correspondant & 
deux surfaces stratigraphiques (S°i, S°„) de la realisation de r6f6rence. 
Un plan (G \ 9 Q m ) de ia description maillSe (G ) representant une des deux surfaces 
stratigraphiques (S k i» $\) de la realisation; (R k ) est obtenu en appliquant aux points 
des plans (G°i. G° n ) de la description de reference les deplacements permettant de 
15 passer de la surface stratigraphique (SV S° n ) de r6f6rence a la surface stratigraphique 
(S k i, S k n ) dans la realisation (R k ). Les autres plans de la description sont obtenus 
ensuite par interpolation. 

Le procede permet de preserver les homologies de structure comme de 
representation entre le$ difftrentes realisations. II evite de recalculer ex nihilo une 

20 description maillee de la reference. j 
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Foumiture de la realisation de reference R°, de la grille 
de reference G° et de la realisation R k 



Choix d'un pilier P° ft h de la 
grille de reference ■ 



Determination dc Vintersection du; pilier avec 
les surfaces de tpit S°i et de base S° n 



8^ 



Caloul du deplaccmcnt dcs intersections lors du passage 
de la realisation de reference R° a la realisation R k 



16 



Oui 
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14' 



Choix d'une surface 
uiterm£diairc S k / 



Choix d'un pilier P k 5 . h 
de la grille 



" Detetrnination de Intersection du ( pilier avec 
la surface intermediaire S k j 




Surface \ ; Oui 
suivante ? 



Fig. 3 
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_ Choix d'un pilicr P k Bfh \ 
de la description maillee 



I 
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Cnleul sur le pilier P° g h dc la description maillee de 
reference de la position des plans G°J de la couche 



Cakul des points des plans G k , sur le pilier P 1 ^, suivant 
les m£mes proportions que sur le pilier P 0 K .h 




Retour a I'etape 32 figure 3 



Fig. 5 
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Choix d'une hypothese 
geologique I 



Determination du pilier E\u le plus 
long de la couche 



54 



Calcul de la position des plans G° r . de la couche sur le pilier 
F'° Bj h correspondant dc la description jrnaillee de reference 



Calcul des points des plans G\ sur le pilier P 1 ^ suivant 
les memes proportions que sur le pilier P 0 ^ 



56 



58 



CalcuJ des points des plans G\ sur les autres pilier 
P k *fc en respectant I'hypothese de parallelisme 



Retour k Tetape 32 figure 3 



Fig. 6 



R.\BfOVClsM0l00M91HOflg<ydoc * 2 1 /US/2002 



8/8 



74 



76 



72 



Caicul des points reguliers 
dc la description mail lee G k 



Choix d'un point irregulicr N° B ,b,i <fc la 
description maillee de reference 



Determination des points reguliers voisins et 
de leur deplacemcnt 



Calc ul du displacement de ia position du point irregulier a 
partir des deplacements des points reguliers voisins 




Fig. 9 



